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In this work I present a numerical study of the Finite Size Scaling (FSS) 
of a correlation length in the framework of the CP^~^ model by means of 
the 1/N expansion. This study has been thought as propedeutical to the 
application of FSS to the measure on the lattice of a new coupling constant 
fx{l/R), defined in terms or rectangular Wilson Loops. I give also a pertur¬ 
bative expansion of fx{l/R) in powers of the corresponding coupling con¬ 
stant in the MS scheme together with some preliminary numerical results 
obtained from the Polyakov ratio and I point out the conceptual problems 
that limit this approach. 


I. INTRODUCTION 

A very important goal of a non perturbative approach to an asymptotically free theory 
such as the lattice is the determination of A;^ in a physical mass unit; this requires a 
very precise measurement of the ’running coupling’ for very large fi. In order 

to perform the measure by a Monte Carlo simulation on a lattice it is useful to work 
in the framework of a new renormalization scheme, with its own renormalized coupling 
constant expandable in powers of ajjg{fi). This can be done dehning a renormalized 
coupling constant in terms of rectangular Wilson loops. It is possible then to use Finite 
Size Scaling techniques to reach very small distances without very large lattices |^. 

In this work I want to test this program in the framework of an abelian gauge theory 
in two dimensions: the CP^~^ model. The choice of this specihc model is justihed 
by the fact that it is a renormalizable gauge theory and, above all, it shows both the 
asymptotic freedom and a conhning potential between a particle-antiparticle pair, just as 
QCD does. Furthermore, CP^~^ is simpler than QCD, being an abelian gauge theory 
in two dimensions and having the possibility of an expansion in powers of 1/N, that can 
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be efficiently used to test the new renormalization scheme and the new approach to the 
problem of computing A;^. 

This paper is organized as follows: in section |1| I introduce the CP^~^ model both in 
the continuum and on the lattice. In section HI apply the FSS technique to the study of 
the correlation length, trying to measure the A-parameter in two different renormalization 
schemes. In section I adopt the dehnition of the running coupling given in |]l| , whose 
perturbative expansion in powers of fjjsifi) is obtained. In section [V] I give some analytical 
and numerical results for the new coupling constant, obtained by means of the Polyakov 
and Creutz Ratio. In section 01 conclude with some considerations about the strategy 
outlined in the present work and I try to single out the ulterior difficulties that come out 
when one wants to apply the same technique to the study of the Creutz Ratio. 


II. THE CP^-^ MODEL 


The CP^ ^ model is a generalization of the non linear cr-models. The bare lagrangian 
of the continuum theory is: 


L [z (x ), df,z (x ), Af, (a:)] 


N 


Dfj,z {x)D^z (x) 


( 1 ) 


where z is & complex A^-vector constrained by the condition z{x)z{x) = 1 and the operator 
is dehned as [| . 

As it is suggested by the mass dimension of the coupling constant /, the theory is 
renormalizable and it shows dimensional trasmutation [^|-Q. 

Furthermore, the theory is expandable in powers of 1/A^; in this framework it is 
possible to show that the massless particles of the ‘standard’ perturbation theory acquire 
mass at the leading order in 1/A^ and they interact by a linear confining potential 

A very good reference for the theory of spin models on the lattice is |P], where we hnd 
a convenient lattice action of the CP^~^ model, the gap equation, the inverse propagators 
and Ag ^ of the effective helds and the dehnition of a correlation length that I 
report here: 
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4sin2| 


Gp (0, 0; L) 

Gp (0,1; L) 


( 2 ) 


^It can be considered as a covariant derivative associated to a U{1) local gauge invariance. 
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Where Gp is the two-point correlation function of gauge invariant operator and Gp is 
its Fourier Transform. We can write Gp{ki, k 2 ] L) as a function of the propagator of the 
effective field a in the framework of the 1/N expansion on a finite lattice: 

Gp (k; L) = iAj -4 (i; L)+0 (i) (3) 

where {ki,k 2 ;L) is given in |^. Using the gap equation on a finite lattice 0 I can 
write the function (/3) in a parametric form (with mp used as a parameter): 


f3 {mp) 


4 sin 2 i 




A-i(0,0;L;mi) 


L2 



(4) 


Let us consider now the inhnite lattice limit of the correlation length so dehned. Ma¬ 
nipulating the formulae in [Q, if is easy to show that ^oc{P) has the following asymptotic 
form: 


^oo (/3) ~ 


^ ^27r/3 

sW 


(5) 


III. FINITE SIZE SCALING OF THE CORRELATION LENGTH 


A complete theory of the Finite Size Scaling is contained in HJTffl, while an explanation 


of the strategy used in what follows and the definition of the various FSS functions can 
be found in |lTl|,p!^. The FSS says that: 


or equivalently 


iP) ^ fp 



^L2 iP) 
a, iP) 


~ (Tg 




( 6 ) 

(7) 


I have used the expressions reported in in the framework of the 1/N expansion, instead 
of a standard Monte Carlo simulations, in order to evaluate numerically the correlation 
length. The numerical computations have been performed by means of a FORTRAN 
code, choosing lattices of sizes varying from L = 20 to L = 2000. 
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FIG. 1. Function for some values of L. 


A. Computation of ahiAT 


The results of these computations are showed on the hgure ^ From these data I can 
extract a measure of the A-parameter in the lattice regularization scheme, dehned by the 
perturbative two-loop solution [| of the renormalization group equation for the coupling 
constant in the large N limit as reported in (j^) 


a (^-^LAT 


( 8 ) 


The difference between aAiAT extracted from the lattice and the theoretical value 8-\/3 
(cfr. d^)) is below 3.2%. An equivalent way to test the asymptotic scaling is to define an 
effective A-parameter | : 


1 




(9) 


^Remember that / = 1/2/3. 

^The function CooiP) is exactly known from |p 
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and to plot the function 


^LAT _ ^ 

^eff 

versus ^oo(/3), as showed in the figure ||. The deviation from the asymptotic behaviour 



FIG. 2. Asymptotic scaling test. 
(H) is less than 2.6 x 10“® for > 90. 

B. Computation of a^{2,ZL) 


and the hgure ^ shows that it has these two asymptotic behaviours: 

lim (Tf (2, zx,) = 1 lim (Tf (2, zx) = 2 (10) 

ZL^oo ^ ^ ' 

that reflect these simple physical considerations: 

1. In the limit zl ^ oo the actual size L of the box is not very important because the 
correlation length ^l{P) of fhe interaction is very small, compared to L; we are very 
close to the inhnite volume limit. 


I have computed the function ae (2, zi) n , according to the scheme outlined in 12 


4 


ZL = 


L 
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2. In the limit Zl —*• 0, on the contrary, the correlation length is rnuch larger 

than L, that becomes the actual scale of the interaction: the observables are very 
sensible to what happens on the border of the box; that is, ^l(/ 3), having the same 
dimension as L, is directly proportional to it. 

If the FSS is correct, the function (J^) should depend only on the ratio and not on 
and L separately. I can check if this is the case by verifying that the several curves 
obtained on the different pairs of lattices approximately superimpose. The results are 
showed on the figure |^. The biggest relative violation of the FSS is about 5.4 x 10“^ for 



Zl=L/4l 

FIG. 3. FSS function crg(2, zl) for some values of L. 


zl = 4.635... 


C. Reconstruction of f^{zL) 


Now, in order to go on with the program outlined in |I^, I need to know the function 
(jg (2, Zl) in every point of a sufficiently large interval; I have two possibilities: 


1. interpolate the points I have obtained numerically using a polynomial of n-th degree, 
with n to be chosen. 
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2. find a suitable function that fits the points at our disposal very well. The theory 
tells us that ct^(2, zl) —> 1 exponentially fast as zl —> oo. Then I can try to use 
a £t function of the form 

N 

{2,zl) = 1 + X! «^exp {-uzl/s) (11) 

n=l 

with the constraint that lim^^^o (2, zl) = 2. 

There are several systematic effects intrinsic in both procedures; I have made some 
checks in order to choose the values of the parameters that minimize them. Anyway, all 
‘spurious’ systematic errors can be made smaller than the FSS violation, which in turn 
constitutes the real physical limit of our procedures. The hgure ^ shows the results of the 
calculation. 
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FIG. 4. FSS function computed in three different ways: 

b.l5—paramters fit of cr^. c.Direct calculation of as . 


a.Interpolation of 


is a scale factor to be chosen. An analysis of the tipical mass unit of the theory with N = oo 
suggests s = 2^/6 as the best choice. 
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D. Perturbative expansion of and 


A perturbative calculation on a finite lattice, along the line sketched in 
following result for the FSS function f^\ 


14] gives the 


U {zl) = R (l + O (4)) (12) 

with R = = 8v^ and d = 0.76077028681... In the large N limit; 


u (zl) = (l + O (4)) (13) 

and 

, , / ln2 9 In^ 2 ^ \ , 

apert (2, Zl) = 2 1^1 - -^Zl - + o {zj)jj (14) 

I checked the range of validity of this perturbative result and I found good agreement for 
all Zl < 1. 

We can do this perturbative test on the function directly. If we look at the shape 
of the function we realize that the number of points in which it is necessary to know 
itself in order to reconstruct is a decreasing function of zl] this implies that a 
systematic error on the numerical computation of will propagate more and more when 
Zl becomes smaller and smaller. This is exactly what I find if I compare the numerically 
reconstructed functions and with the exact perturbative result ([I3|). It is worth 
noticing that I can explain this effect quantitatively quite well invoking the FSS violation 
studied at the end of the paragraph |III B| . 


E. An asymptotic scaling test 

It is worth noticing that in the large N limit the FSS function f^{zL) is exactly known 


in terms of the two functions 15 


Fi{z) = -log TT H-— + 7 e+ 

z 2 


+2vrE 


+ 


2vr E 


\ + 2)2 2Tm^ 

1 1 
V^vr%2"+z2 exp + z'^'^ - 1 


(15) 






















( 16 ) 


F2{^)= E 


coth 


V47r^n^+2:^ 

2 


^=_oo 1 “ 4n^ \/ATT^n?~+^ 

so we can obtain the fnnction f^izi) in a parametric form: 

^ TZil 


(17) 


As already seen at the beginning of this section, the asymptotic behaviour of this expres¬ 
sion is given by (^); in order to test the speed of approach to the asymptotic scaling I 


define from equation (|T^) the ‘asymptotic function’: 

= ^FSS 


47r 


zl 


(18) 


where Apss = Re ^ and d = 0.76077028681.... I define also an effective Ai? 5 ' 5 -parameter: 

47r 

^FSsizL) = f^{zL)zLe^ (19) 

and I plot in the figure ^ the function 


^^FSS 

Afss 


- 1 


( 20 ) 


that shows the deviation from the asymptotic scaling; this is less than 3.5 x 10 ^ for 
zl < 0.036. 


IV. THE ‘RUNNING COUPLING’ FROM PERTURBATION THEORY 

Following the procedure described in |II], let us consider a rectangular Wilson loop 
{R X T) 


Wr(R,T) = 


( 21 ) 


It is possible to give a definition of an interaction force Xt{R) in which a partial derivative 
with respect to T substitutes the lim^^oo of the standard dehnition |1l6tp!7| : 


x{R,T) = 


\nW{R,T) 
dRdT ■ 


( 22 ) 
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FIG. 5. Asymptotic scaling test for the FSS function 


The perturbative expansion, in the continuum theory, of the interaction force Xt{R) so 
defined is: 

x¥(R) = {l + 2''o !>■>(/*«) + 0(^)1 + 0(fe(f‘)")} (23) 

where 


(2vr)^ qi go 


(24) 
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and 


In a: 


1 

TT 

+ 


h{x) = 


7s 


1 + 




-B(x) H —- arctan x In x + 


x^ 


1 + x^ 


/•i (in 
/o M 


■ arctan 


MX 


x^ 


1 + x2(l — u) ^ 


^ + i|.i(l + x^)(^ + i)] + 

+ arctan (--) \ + 


1 + x^ 


u 


with 


and 


X 


+ — / dtt^\n(^ — 1 

TT JO \t^ 


+ 


(l + x^t^)^ (x^ + t^)^ 


B{x) = g{x) + ^g{-) 


x^ 


X 


9{x) 


( 7 £; - 1) arctan- 1e ——^ 

X 1 + x^ 


1 , X 

arctan — In-- 

X 1 + x"* 


( 28 ) 


(29) 


(30) 


These fnnctions have the following x —>■ oo limit: 

lim /i(x) = —- lim / 2 (x) = — - - (31) 

Now it is possible to compnte the fnnction (|^) nnmerically. The resnlts are reported in 
the table |. 


TABLE I. Numerical value of the function a(x) for some values of x. 


X 

a{x) 

X 

a{x) 

1 

-0.289942... 

4 

-0.091605... 

5/4 

-0.450157... 

5 

0.007810... 

4/3 

-0.472398... 

10 

0.238997... 

3/2 

-0.488345... 

50 

0.478908... 

5/3 

-0.479687... 

100 

0.520289... 

7/4 

-0.469696... 

200 

0.544656... 

2 

-0.427594... 

300 

0.553868... 

3 

-0.234557... 

oo 

0.577216... 


I can now dehne a rnnning conpling constant 


x{R,T) 


N-1 

271 


/l(x) 


txii/Rf 



(32) 
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and a large N-rescaled coupling constant f^iX/R) = Nt^iX/R)/2 


,Ul/RY 

x{R,T) = -^7:^h{x)- 


( 33 ) 


71 ' ' IP 

The definition and the perturbative result (^) allow to establish a link between 
4 (T) and 


/ 1 \ Z 

*“’ ( (r) “ ('*)d ^ 1“(^) + '*■^1 ® (*MS (pf) } 


(34) 


where a(x) determines the connection between the A-parameters in the two renormaliza¬ 
tion schemes: 

a(x) = In (35) 

^Ws 


V. ANALYTICAL AND NUMERICAL RESULTS FROM THE POLYAKOV 

RATIO 


The Wilson line or Polyakov loop is strictly related to the static quark-antiquark 
potential [0. I can take the interaction potential in the continuum theory for the case 
X —>■ oo from [01: 


(ttt 2rr 2 

V{R) = -^mlR + cos{kR) 


^ ,,, 127rmn 


dk 


(36) 


where 




2ir 




Y(k) = xll 


ArriQ 

1^ 


(37) 


After having rotated the integration contour in the complex plane and other algebraic 
manipulations I obtain 


V(m - , 2ci 27r p 

V[,R)-—mQR+—mQ + — e 

iV iV I\ J 2mo I 


Y 


(7rY)2 + (Vln^ -2)^ 


from which I can extract the interaction force 


TT 


NF(R)^ = TT H— / dte 
3 J2 


tY 


{ttYY + (Vlni 


i+y _ 9 

-y ^ 


-R 

2 = <h(r = -) 


(38) 


( 39 ) 


where in (|38|) and (|39D Y{t) = Jl — 
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The formula (|3^ can be used to compute numerically the reference continuum quantity 
to be compared with the lattice results; it can also be expanded in the regime r —0 
(perturbative or scaling region) obtaining the scaling behaviour of the interaction force: 


$(r) = 


TT 


t^Ie 


+ 0 {- 


2 j -2 ^2 ^ • 

r r r 

This result can be rewritten with the same precision as 


<h(r) = 


TT 


2 In^ 

e^Er 


l + 0 ( 


J .2 vT ' 
r 


( 40 ) 


( 41 ) 


It is possible to extract from (^) the running coupling constant in the scheme outlined 
in the previous paragraph in the limit x ^ oo 


focir) = 


71 


In 


vT 

e^Er 


( 42 ) 


This expression is in perfect agreement with the prediction of the 1-loop perturbative 
renormalization group 


<h(r) 


TT 1 
2 In^ 


1 

AooK 


( 43 ) 


if the A-parameter is chosen consistently with the renormalization scheme introduced 
before Q : 


Aoo-R = = e'^^r/VQ 


(44) 


The difference between the exact continuum curve and the asymptotic one is below 
5% for r < 0.25. One could then measure the interaction force on a finite lattice by means 
of the so-called Polyakov Ratio, for several values of < 0.25, keeping zl fixed, and then 
fit the values with equation (|43|) in order to extract the A-parameter. It is so possible to 
understand how fast the measures on the lattice approach the continuum limit (|^) when 
the lattice spacing goes to 0 . 

The Polyakov Ratio on a finite lattice with periodic boundary conditions is given by 


1 W{R,L) _ 1 l^sin 

^ L ’^1T(R-1,L) 2NL^ 


ki[2R—l)7r 


fclTT 


-ki A^g){ki,0) 


(45) 

k ^ 

in the leading order in the 1/N expansion, where A( 0 )(fc) is the lattice propagator of the 
gauge held 9 reported in [^. 


®See (Ht) 
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FIG. 6. Systematic effect oc a/R on the value of the Polyakov ratio on a finite lattices. The 


continuum line is the continuum interaction force 40 


A numerical calculation of the Polyakov Ratio has been performed on two lattices 
{L = 1000 and L = 2000) and for two keeping zl = 19.455... hxed. The results are 
showed on the hgure H, where we can see the presence of a systematic effect oc a/i? for 
r = R/^l < 0.25; when r > 0.25 other effects due to periodicity start to play a much 
more important role P,|T^ . 


VI. CONCLUSIONS AND OUTLOOK 


The study of the correlation length has clarihed how to use the FSS technique to 
estrapolate the hnite volume measures to inhnite volume results. In this framework it 
could be interesting to apply the same technique to a basic observable: the Creutz Ratio 
derived from rectangular Wilson loops. Even if this seems to be just a mere step-by-step 
execution of the program already outlined in [M for the correlation length, there are 
some points that should be underlined: 


1. the presence of the additional scale R implies that the FSS function depends on two 
variables, instead of only one. 

2. a systematic effect oc (analogous to that found in the case of the Polyakov 

Ratio) is present in the Wilson loop case, too. 
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3. in the special case of an abelian lattice gauge theory dehned on a torus, as our 
version of the lattice CP^~^ model is, there are some kinematical effects that have 


to be considered in order to obtain significant results |l9 


4. the special shape of the FSS function cr^ (2, zl) imposes another limitation on the 
region where the FSS technique can be applied efficiently: let us define 


t(/ 3 , R, X, L) = x(P, R, xR, L) 


( 46 ) 


$ obeys to a FSS law similar to (j^) seen in the section mi for the correlation length: 

( 47 ) 


* (i, /?. R, L) = 9J, h. 


The counterpart of the function is 

Sa,b{zL,rL) = lim 


<h(/9, bR, X, aL) 
d^oo <k(/3, -R, X, L) 


(48) 


where 


L _ R 

JM 


( 49 ) 


The limit (5 ^ oo means simply that we are in the scaling region. If the FSS law 
(0) is valid and choosing a = 2, 6 = 1 


S2,l{.ZL,rL) 


g^{x,z'L,r'y^) 

g^ix,ZL,rL) 


<F(/3, R, X, 2L) 

Hi^,R,x,L) 


( 50 ) 


with 


^zl , _ Tl 

^ a^{2,ZL) 


(51) 


Once I have measured the function S 2 ,i(^l, ?’l) in a large enough region of the plane 
(zLyri), I can try to reconstruct the function g<s>{x, Zl^tl) by: 


n 

g^{x,ZL,rL) ^ n - 71irRJK3<s>ix,ZL^ 

j=l S2,i(4 ) 



( 52 ) 


with 


15 















(53) 


ji) _ ‘‘n ji) _ n 

A severe limitation to this computational scheme is represented by the practical 
impossibility to reach very small values of From hgure |] we see that 

(jg ~ 2 for L/^l{(3) < 1. When we try to compute for a L/^l{(3) quite small, 

in (^) stays very close to 2 even for J very large; now, for j < J, while 
tl is reduced to a half at each step, zl stays about constant (remember (|^ ). I 
can do two things, in order to lessen I can lessen R, and, once I have reached 
the lowest limit for R, since a is hnite, I must increase and this means that I 

have to increase L too, to keep zl constant. In conclusion, if I want to compute the 
function gq,{x, zl^tl) for a very small zl I need to simulate on very large lattices, 
or, in other words, a upper limit on the size L of the lattice implies a lower limit on 
Zl, as far as the measure of g^{x, zl^ti) concerns. 

These considerations leads to the conclusion that the application of the FSS technique 
to the Wilson loop needs a careful study of all the systematic effects. Furthermore, the 
problem outlined in the last point above seems to limit the efficiency of the method in 
the region where also the traditional techniques fail. The crucial point is represented by 
the exact form of ct^(2, zl) when zl —^ 0. Once I have hxed the largest dimension Lmax 
of a usable lattice, the lowest limit at which it is possible to compute the function 

g<s>{x, Zl, tl) is different from 0 and it depends on the specihc dehnition of the correlation 
length chosen. The only thing we can do is to hnd a correlation length that minimizes 

(min) 
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